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A  parabolic  equation  model  for  scattering  in  the  ocean 
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The  small-angle-of-propagation  limit  and  the  method  of  matched  asymptotics  are  applied  to 
derive  an  efficient  model  for  solving  realistic  underwater  acoustics  problems  involving  both 
propagation  and  scattering  from  a  submerged  object.  The  propagation  and  scattering  aspects  of 
the  waveguide  scattering  problem  are  decoupled  by  approximating  the  waveguide  Green’s 
function  on  the  surface  of  the  scatterer.  For  low  frequencies,  the  small-angle  limit  also  allows 
one  to  approximate  the  incident  field  with  a  horizontally  propagating  plane  wave  and  the 
scattered  field  with  an  azimuthally  specular  point-source  field.  With  these  approximations, 
scattering  calculations  can  be  performed  efficiently  in  the  time  domain.  Calculations  involving 
the  three-dimensional  parabolic  equation  and  the  time-domain  parabolic  equation  are 
presented. 

PACS  numbers:  43.30.Gv,  43. 30. Bp,  43. 30. Dr 


INTRODUCTION 

Several  models  have  been  developed  for  solving  prob¬ 
lems  in  both  underwater  propagation1  and  free-space  scat¬ 
terings  Yet  these  remain  highly  active  areas  of  research. 
Thus  it  is  not  surprising  that  model  development  for  scatter¬ 
ing  in  a  waveguide,  a  problem  involving  both  propagation 
and  scattering,  has  only  recently  begun.  The  few  models  that 
have  been  developed  ’  do  not  take  advantage  of  asymptotic 
limits.  Thus  these  models  lack  both  simplicity  and  efficiency 
and  are  difficult  to  apply  to  realistic  problems. 

The  small-angle  limit  (nearly  horizontal  rays)  is  per¬ 
haps  the  most  useful  asymptotic  limit  in  underwater  acous¬ 
tics.  This  limit  is  used  to  derive  the  parabolic  equation 
(PE),'’  which  gives  accurate  and  efficient  results  for  propa¬ 
gation  problems  in  complicated  ocean  environments.  The 
three-dimensional  PE  (3DPEF  ’’  is  required  in  general. 
However,  the  two-dimensional  PE  (2DPE)  is  valid  when¬ 
ever  azimuthal  diffraction  is  negligible. The  inverse  Four¬ 
ier  transform  is  applied  to  the  PE  to  derive  the  time-domain 
PE  (TDPE)."  "’  which  is  an  accurate  and  efficient  pulse 
propagation  model.  If  the  source  and  receiver  are  sufficient¬ 
ly  far  from  the  scatterer,  these  models  can  be  used  to  propa¬ 
gate  the  incident  field  to  the  scatterer  and  the  scattered  field 
to  the  receiver. 

For  a  point  source,  the  nearfield  effects  of  attenuation 
and  variations  in  sound  speed  and  density  do  not  affect  the 
farfield.17  Since  a  scatterer  behaves  like  a  source,  this  sug¬ 
gests  that  the  far  waveguide  scattered  field  is  not  affected  by 
refraction,  reflection,  and  dissipation  near  the  scatterer.  We 
demonstrate  this  using  the  results  of  Ref.  17  and  the  method 
of  matched  asymptotics."4  We  show  that  the  scattering  and 
propagation  aspects  of  the  waveguide  scattering  problem  de¬ 
couple  into  the  homogeneous  half-space  scattering  problem 
and  simple  propagation  problems. 

The  half-space  scattering  problem  is  further  simplified 
foi  low  frequencies  using  the  small-angle  limit.  The  incident 
field  ean  be  approximated  by  a  plane  wave,  and  the  scattered 
field  can  be  approximated  by  the  field  due  to  an  azimuthally 
specular  point  source.  These  approximations  allow  efficient 


pulse  scattering  calculations.  Scattering  calculations  involv¬ 
ing  the  TDPE  and  the  3DPE  are  presented. 

1.  THE  FARFIELD  APPROXIMATION 

Cylindrical  coordinates  are  used  with  r  being  the  depth 
below  the  ocean  surface.  0  being  the  azimuth  angle,  and  r 
being  the  horizontal  distance  from  the  center  of  the  scatterer 
T  of  length  L  at  z  =  z„.  The  point  source  Si  is  located  at 
( r.z.O )  =  (r,^,M,).  The  coordinate  system  and  geometry 
are  illustrated  in  Fig.  1 .  The  ocean  depth  is  denoted  by  d.  The 
acoustic  pressure  p  satisfies  the  pressure  release  boundary 
condition  at  the  ocean  surface  and  the  outgoing  radiation 
condition  at  infinity.  For  the  time-harmonic  problem  with 
circular  frequency  to,  the  complex  pressure  P  is  defined  by 

p(r,z,0,t)  =  P(r.z,0)c\ p(  -  icJt).  ( 1 ) 

where  l  is  time.  Away  from  T  and  S l.  P  is  assumed  to  satisfy 
the  reduced  wave  equation  for  a  medium  with  variable  den¬ 
sity"' 

/jV%\ ( \/p)VP\  4-  A'  P  =  0.  (2) 

where  p  is  the  density  and  the  complex  wavenumber 
K  =  k  +  ii/ft  k  is  used  to  model  attenuation  that  depends 
linearly  on  frequency. The  real  wavenumber  is  k  =  oj/c, 
where  c  is  the  sound  speed.  The  reference  wavenumber  is 
ku  =  k  (z„).  the  reference  sound  speed  is  o„  =  £•(;,,).  the  sedi¬ 
ment  attenuation  in  decibels  per  wavelength  ( dB/x  )  is  ft. 
and  »/  =  (40-  log,,,  e) 

Rays  that  propagate  within  the  critical  angle  6,  from 
horizontal  are  called  trapped  rays.  We  define  e  =  tan’  6, 
and  assume  that  1 .  The  critical  range  r  ,  which  is  defined 
by  r  tan  cA,  =  d  —  z„,  is  the  range  at  which  the  steepest 
trapped  ray  from  the  center  of  F  intersects  the  ocean  bottom. 
To  justify  using  the  PE  approximation,  we  assume  that  range 
and  azimuth  variations  in  the  environment  are  weak.  The 
incident  field  P,  is  defined  to  be  the  field  due  to  Si  in  the 
absence  of  F,  and  the  scattered  field  P  is  defined  by 
P  —  P,  4  P  .  The  waveguide  Green’s  function  6’(x,x’ )  satis¬ 
fies  the  pressure  release  boundary  condition  at  the  ocean 
surface,  the  outgoing  radiation  condition  at  infinity  ,  and 
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FIG.  I.  Coordinate  systems  and  geometry:  (a)  View  from  the  side;  (b)  view 
from  above  the  ocean  surface. 


pV-l(l//3)VGl  +  K'G=  -4  jrb(x-x’),  (3) 

x  =  (x,y,z)  =  (/-cos  (Arsin  6,z).  (4) 

Our  definition  of  G  differs  from  the  usual  definition  of  the 
Green's  function  by  a  factor  of  4^-.  The  scattered  field  has  the 
representation 

dP  \ 

- (x')G(x.x')  \dA  '.  (5) 

dn  ) 

where  r?  /dn  is  the  outward  normal  derivative.  G  is  the  field 
d  te  to  a  point  source  at  x  =  x'.  and  dG  /dn  can  be  approxi¬ 
mated  with  arbitrary  accuracy  by  a  pair  of  point  sources  near 
x  —  x'  Thus  P  can  be  approximated  with  arbitrary  accura¬ 
cy  as  the  field  due  to  a  superposition  of  point  sources  on  and 
near  dV. 


In  the  homogeneous  free  space  with  c  =  c„  and  0  =  0, 
the  integral  representation  can  be  used  to  obtain  the  follow¬ 
ing  farfield  expansion  for  the  free-space  scattered  field  ( P , ); 
for  an  incident  plane  wave  of  unit  amplitude21 


(P,)  ,(x)  ~  Dj(  6,<J> )  G,  ( x,x() ) ,  (6) 

G,(x.x')  =  [exp(zA-„]x  -  x'|)]/|x  -  x’j,  (7) 

x„  =  (0,0,z„),  (8) 

<b  =  tan  1 1  (z  -  z„)/rj.  (9) 


For  large  r.  dH  /dr  =  0(k,,e),  dll /dz  =  0(kue'  :),  and 
dH /d04kltr  ( Ref.  8).  For  the  case  k„L  =  0(1),  the  weakly 
varying  coefficient  function  H(x.x')  can  be  replaced  by  the 
constant  H(x,x„)  for  x’  e  d T  to  leading  order  in  e.  With  this 
approximation,  the  integrand  in  Eq.  ( 5 )  is  reduced  to  a  form 
in  which  the  analysis  of  Ref.  2 1  can  be  applied  for  /■>  L.  We 
obtain 

\  ( .V  -  x'r  -f  (y  -  y')-  sc  -  x'-r.  (11) 

G(x,x')  s=G(x.x„)exp(  -  ikt)x '•?),  (12) 

DG 

— (x.x')  sr  -  /A,//*AG(x,x„)exp(  -  ikux'-r),  (13) 

dn 

where  r  =  (cos  O.sin  0,0)  and  n  is  the  outward  normal  unit 
vector  on  r?r.  Substituting  Eqs.  (12)  and  (13)  into  Eq.  (5). 
we  obtain  the  waveguide  analog  of  Eq.  (6)  for  large  r: 

P,(x)-~D(0)G(x,xo),  (14) 

where 


D(0)  =  -  JJ(/A„/Vr/,(x')  +  -y^  (x  )^ 

Of 

Xexp(  —  ik(lx'"r)dA 


(15) 


Thus  the  scattered  field  can  be  approximated  by  a  specular- 
point-source  field  for  r^>L.  This  result,  which  was  obtained 
for  a  rigid  sphere  in  Ref.  5,  is  in  agreement  with  some  obser¬ 
vations  for  multipole  sources.'  Since  Eq.  (14)  is  a  leading- 
order  approximation,  it  should  be  most  accurate  in  the  for¬ 
ward  and  specular  scattering  directions,  where  P  is  most 
intense.  In  other  directions,  P \  is  relatively  w  eak.  and  higher- 
order  asymptotic  terms  may  be  important. 


II.  THE  MATCHED  ASYMPTOTIC  SOLUTION 

The  matched  asymptotic  solution  of  a  differential  equa¬ 
tion  is  composed  of  the  inner  solution,  which  is  valid  in  the 
inner  region  of  the  domain,  and  the  outer  solution,  which  is 
valid  in  the  outer  region  of  the  domain  (sometimes  an  inter¬ 
mediate  region  is  required).  Different  terms  of  the  differen¬ 
tial  equation  are  important  in  the  two  regions.  For  example, 
we  consider  the  matched  asymptotic  solution  of  Eq.  ( 3 )  A 
leading-order  inner  solution  is  obtained  by  assuming  that  the 
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inner  region  near  fl  is  homogeneous  with  c  =  c(t  and  ft  =  0 
(Ref.  17).  In  the  outer  region  far  from  fl  (At„r>l),  the 
asymptotic  expression  for  the  Hankel  functions  appearing  in 
the  normal  mode  solution3'24  of  Eq.  (3)  is  valid.  Thus  Eq. 
(3)  may  be  replaced  with  the  PE  in  the  outer  region. 

If  the  inner  and  outer  regions  overlap,  the  following 
matched  asymptotic  solution  of  Eq.  (3)  is  obtained: 


f  G,n  r  <  r0, 
1  G„„,  r>r,„ 


(16) 


where  r0</-(  (Ref.  17).  The  inner  solution  Gm  is  the  homo¬ 
geneous  half-space  Green's  function,  which  satisfies  the 
pi’ssure  release  boundary  rondition  at  z  =  0,  the  outgorng 
radiation  condition  at  infinity,  and 

V3Gin  +  k- G„,  =  -  4m5(x  -  x'),  (17) 

Gjn  (x>x')  =  Gf(x,x')  -  Gf(x*,x),  ( 18) 

x*  =  (x,y,  —  z).  (19) 

For  many  problems,  G,n  excites  the  propagating  normal 
modes  more  accurately  than  the  Gaussian  PE  starter. 17  Al¬ 
though  G,„  breaks  down  for  k(,(z  —  d)  y  1,  this  does  not  af¬ 
fect  the  accuracy  of  G„ul  because  rays  that  enter  the  sediment 
in  the  inner  region  are  not  trapped.  In  the  outer  region  r  >  r„, 
the  outer  solution  satisfies 


pV*(— VG,„")  +  =°,  (20) 

with  the  boundary  condition  Gtml  —  Gm  at  r  —  r„.  The  PE 
can  be  used  to  solve  Eq.  (20)  if  £-„/•„>  1. 

Treating  the  region  near  a  point  source  as  a  lossless, 
homogeneous  half-space  causes  negligible  errors  in  the  far- 
field  because  refraction  in  the  ocean  is  weak,  attenuation  in 
the  ocean  bottom  is  weak,  and  rays  that  reflect  from  the 
ocean  bottom  near  the  source  propagate  at  large  vertical  an¬ 
gles  and  are  not  trapped  by  the  oceanic  waveguide.  Since  T 
behaves  like  a  collection  of  point  sources,  one  would  expect 
that  an  analogous  matched  asymptotic  solution  that  in¬ 
volves  the  solution  of  the  half-space  scattering  problem  ex¬ 
ists  for  the  waveguide  scattering  problem. 

The  waveguide  scattering  problem  can  be  posed  as  an 
integral  equation  involving  G.  We  do  not  write  down  the 
integral  equation  because  its  form  depends  on  the  properties 
of  T.  An  asymptotic  solution  of  the  integral  equation  is  ob¬ 
tained  by  approximating  G  with  G„,  on  <?T,  which  gives  the 
integral  equation  for  the  homogeneous  half-space  scattering 
problem  having  the  solution  (Ps  ),„  with  the  total  field  Pm. 
This  approach  is  an  alternative  to  the  Kirchhoff  method, 
which  is  used  to  solve  scattering  integral  equations  by  ap¬ 
proximating  Ton  r?r.  Thus  we  obtain  the  following  matched 


asymptotic  solution  for  T, : 

f(^)„.  7<r„, 

it/7,) . 

(21) 

‘  *</*,) . )  F  K  ?(P  ),„„  -0. 

0  ) 

(22) 

with  the  boundary  condition  ( P  ),„„  -  ( T  ),„  at  r  ~  /-„.  The 
PF.  can  be  used  to  solve  Eq.  ( 22 )  if  1 .  Decoupling  the 
waveguide  scattering  problem  is  a  useful  simplification  be¬ 
cause  models  for  solving  the  half-space  scattering  problem 


and  the  propagation  problem  in  the  outer  region  have  been 
developed. 

The  matched  asymptotic  solution  for  a  point  source  is 
very  accurate  in  the  outer  region,17  and  T  behaves  like  a 
collection  of  point  sources.  Unlike  a  point  source,  however, 
T  can  induce  multiple  reflections  between  itself  and  the 
ocean  surface  and  bottom.  The  matched  asymptotic  solution 
accounts  for  multiple  reflections  between  f  and  the  ocean 
surface.  Unless  T  is  very  close  to  the  ocean  bottom,  rays  that 
travel  between  two  points  on  dV  via  bottom  reflection  meet 
the  ocean  bottom  at  nearly  normal  incidence  for  which  the 
reflection  coefficient  is  small.  Furthermore,  many  of  the  rays 
involved  in  multiple  reflections  between  T  and  the  ocean 
bottom  propagate  away  from  F  with  large  angles  and  are  not 
trapped  within  the  water  column.  Thus  the  matched  asymp¬ 
totic  solution  for  T,  should  be  very  accurate  in  the  outer 
region. 

Although  (F,  )„,  breaks  down  for  k„(z  —  d)  >  1,  the  fol¬ 
lowing  asymptotic  solution  is  valid  for  k„{z  —  d)  >  1: 


T,(x) - L  f  fY^(x,x')T,„(x') 

477  J  J  V  dn 

«,x'))c 


<*r 


-(x')G(x,x')  jdA 


(23) 


Since  the  waveguide  Green’s  function  is  used  in  Eq.  (23), 
this  representation  accounts  for  refraction,  reflection,  and 
attenuation  in  the  region  r<  r„.  Since  Eq.  ( 16)  is  very  accur- 
age  in  the  outer  region,  however,  it  is  not  beneficial  to  use  Eq. 
(23)  in  the  outer  region. 

Let  us  outline  the  sequence  of  steps  involved  in  the  solu¬ 
tion  of  scattering  problems  with  the  model.  The  incident 
field  is  computed  with  a  propagation  model  such  as  the  PE. 
The  3DPE  must  be  used  if  azimuthal  diffraction  is  impor¬ 
tant.  A  scattering  model  is  then  used  to  compute  (T,  )„,  at 
r  =  r„.  Since  the  accuracy  of  (/*,),„  improves  as  r„  decreases, 
it  is  desirable  to  choose  ru  as  small  as  possible.  Since  many 
scattering  models  assume  an  incident  plane  wave,  it  may  be 
useful  to  use  the  waveguide  normal  modes  to  decompose  P. 
into  plane  waves  as  discussed  in  Ref.  5.  Although  the  half¬ 
space  scattering  problem  is  much  easier  than  the  waveguide 
scattering  problem,  it  is  more  difficult  than  the  free-space 
scattering  problem.  Multiple  reflections  between  T  and  the 
ocean  surface  give  rise  to  mostly  wide-angle  rays  if  T  is  not 
too  close  to  the  ocean  surface.  Thus  multiple  reflections  may 
be  ignored  if  z(l  is  sufficiently  large.  This  aproximation. 
which  was  used  in  Ref.  5,  is  equivalent  to  imaging  the  free- 
space  scattered  field  to  obtain 

(/*,),„  (x)==(Tj,(x)  -  (T,  )/(x#  ).  (24) 


A  propagation  model  is  used  to  solve  the  outer  problem  for 
r>  rtl.  The  PE  may  be  used  if  kltr„$>  1.  The  3DPE  is  required, 
in  general,  because  P  and  possibly  c,p,  and  ft  depend  on  0. 

In  stratified  environments,  however,  azimuthal  diffraction 
diminishes  rapidly  with  range,  and  the  3DPE  is  required 
only  in  the  intermediate  region  r„  <  r  <  r,.  The  2DPE  is  valid  3 
for  r>r,. 
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III.  LOW-FREQUENCY  APPROXIMATIONS 
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For  the  remaining  discussion,  we  assume  that 
koL  =  0(1).  The  results  in  the  remaining  sections  are  valid 
for  scatterers  of  arbitrary  composition,  shape,  location,  and 
orientation.  For  the  examples,  however,  T  is  taken  to  be  a 
rigid  spheroid  of  length  L  and  girth  ttB,  where  k(tB  =  0(1), 
with  its  major  axis  parallel  to  the  x  axis.  Since  fl  is  far  from  T 
and  range  and  azimuth  variations  are  weak,  the  PE  can  be 
used  to  compute 

P,  =  f/,(r',z,0')exp(;A:(/),  (25) 

where  r  is  the  horizontal  distance  from  11  and  6 '  is  the  azi¬ 
muth  angle  about  11.  Far  from  11,  dU,/dr' 
=  O(k0e),  dUl/dz  =  0(k„e'12),  and  dV,/dO'^k,/  (Ref. 
8 ) .  Thus  P,  can  be  approximated  to  leading  order  on  dY  by  a 
plane  wave 

P,  ~  U,  (r,,zo,0,  +  tt) 

Xexp[/A:0(r,  —  x  cos  6,  —  y  sin  0, )  ] .  (26) 

This  is  a  useful  approximation  because  many  scattering 
models  assume  an  incident  plane  wave.  If  B4L ,  the  accura¬ 
cy  of  Eq.  (26)  depends  on  the  orientation  of  T.  Higher  fre¬ 
quencies  are  permitted  for  broadside  incidence  than  for  end- 
on  incidence. 

A  standard  approach  for  matching  the  inner  and  outer 
solutions  of  a  matched  asymptotic  solution  is  to  consider  the 
outer  expansion  of  the  inner  solution.25  For  trapped  rays, 
D/(0,d>)  ~Dj  (OrO )  in  the  limit  e  — 0.  Thus  we  use  Eq.  (6) 
with  tp  =  0  and  Eqs.  (18)  and  (24)  to  obtain  the  outer  ex¬ 
pansion  of  (P,  ),„: 

(P,  C  =  (/,  (r,,zo,0,  +  rr ) exp ( ikt)rt )  /),( 0,0) 

X  Gm  (x,x„).  (27) 

For  the  case  r,  =  rn  (no  intermediate  region  required),  we 
use  the  fact  that  G"“'  =  G,„  and  the  matching  principle25  to 
obtain  the  specular-point-source  approximation 

(P,)„„,  =  P,(r,,zo,0,  +  v)Df(d,0)G(x,x{i).  (28) 

This  approximation  produces  a  significant  gain  in  efficiency 
because  Df,  which  is  difficult  to  compute,  appears  in  Eq. 
(28)  only  for  <b  =  0.  We  deduce  that  D(0)  oc  Df(d, 0). 

To  investigate  the  accuracy  of  Eq.  ( 28 ) ,  we  consider  the 
problem  described  in  Table  I,  where  subscripts  w  and  b  de¬ 
note  water  and  bottom  values.  The  free-space  radiation  pat¬ 
tern  20  log,,, 10,(0, 0)|  appears  in  Fig.  2  for  large  r  and  for 
r  =  rt .  Since  the  radiation  pattern  has  nearly  settled  down  at 
Eq-  (28)  should  be  a  good  approximation  for  this 
problem.  Since  rlt  is  substantially  less  than  r,.,  Eq.  (16)  is 
very  accurate.  However,  it  is  apparent  from  Fig.  3  that  the 


TABU'-  I  Data  for  the  frequency-domain  problem. 


ot 

200  ITS  ' 

d 

225  m 

-  15fX)  m/s 

<:f, 

!MX)  m/s 

/’/ 

1,5  g/L-m' 

II  „ 

0.5  dH//i 

1. 

50  m 

R 

10  m 

(b, 

--  20.4’ 

r, 

-  202.1  m 

r,  - 

l(XX)  m 

rtt 

-  75  m 

z„ 

--  150  m 

r 

5  km 

z , 

50  m 

0, 

-  135“ 

z,  - 

50  m 

FIG.  2.  The  free-space  radiation  pattern  at  r  —  r,  (solid  curve)  and  for 
large  r  (dashed  curve). 


free-space  radiation  pattern  at  r  =  r„  has  not  settled  down. 
Thus  we  use  the  3DPE  for  /■„</•<  r, . 

We  used  Eq.  (26)  to  simplify  P,  to  a  plane  wave  and  Eq. 
(24)  and  the  T-matrix  method2'’-27  to  compute  D,  and 
(P,  ),„.  Transmission  loss  for  P,  and  P,  computed  with  Eq. 
(28)  and  with  the  PE  using  (P,  )„,  as  an  initial  condition 
appears  in  Fig.  4.  From  the  data  in  Table  II,  we  observe  that 
the  approximation  Df(6,d>)^Dl(0, 0)  is  not  accurate  for 
<t>  <d>,  in  the  directions  0  =  0°,  180°  and  225°.  This  explains 
the  differences  in  the  solutions  for  P,  appearing  in  Fig.  4  for 
these  directions.  The  solutions  are  in  excellent  agreement  in 
the  other  directions  including  the  forward  and  specular  di¬ 
rections. 
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Loss  (<IH  re  tin)  Loss  <d II  re  lin)  Loss  (dfi  r»“  Iml  Loss  (dB  re  lmj 


TABLE  II.  Behavior  of  20  log,„|Dy  ( 0,4>)\ 


o 

II 

II 

11 

q 

i>  =  15’ 

d>  =  20° 

0  = 

0“ 

11.69 

12.14 

13.40 

15.27 

17.50 

6  = 

45° 

25.68 

25.70 

25.75 

25.80 

25.78 

e  = 

90° 

13.12 

13.08 

12.96 

12.76 

12.48 

0  = 

135° 

5.37 

5.35 

5.25 

4.95 

4.23 

0  = 

180° 

-  2.07 

-  2.27 

-  2.93 

-  4.32 

-  6.72 

0  = 

225° 

-  8.15 

-  8.68 

-  10.49 

-  14.67 

-  20.07 

0  = 

270° 

-  0.06 

-  0.04 

0.02 

0.10 

0  20 

0  = 

315° 

26. 1 5 

26.03 

25.68 

25  11 

24.35 

where  E„  =  D  if  Eq.  (14)  holds.  The  normal  modes  i /<„  and 
eigenvalues  A„  are  defined  in  Refs.  23  and  24.  We  define  the 
modal  radiation  pattern  associated  with  if',,  to  be 
20  log„,|E„  (0)|.  For  the  stratified  environment  described  in 
Table  I,  modes  I,  3,  6,  and  10  propagate  with  angles  of  ap¬ 
proximately  1.82°.  5.48°,  11.04°,  and  18.62°.  The  radiation 
patterns  for  these  modes  at  r  =  r,  appear  in  Fig.  5  with  the 
far  free-space  radiation  pattern.  The  agreement  with  the 
free-space  pattern,  which  is  excellent  for  mode  1.  decreases 
with  mode  number.  This  is  consistent  with  the  small-angle 
asymptotics. 


The  following  asymptotic  expression  is  valid  in  a  strati¬ 
fied  environment  for  A0r>  1: 

/\  ~  rri^  (z(l)il>„  (z)H o"(A„r),  (29) 


IV.  LOW-FREQUENCY  PULSE  SCATTERING 

To  derive  a  matched  asymptotic  solution  for  pulsed 
sources,  the  Fourier  transform  is  inverted  in  Eq.  (27).  The 
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FIG.  5.  The  far  free-space  radiation  pal  tern  ( dashed  eurves )  and  I  lie  modal  radial  ion  pal  terns  I  solid  curves )  lor  (  a  I  mode  I .  (  h  I  mode  .V  (el  mode  (s.  and 
(d)  mode  10 
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TABLE  HI.  Data  for  the  time-domain  problem. 


<0  =  50  ir  s  ' 

d  =  250  m 

clt 

=  1 500  m/s 

c„  =  1600  m/s 

Pi,  =  1.5  g/cm' 

p, 

=  0.5  dB//i 

L  =  50  m 

B=  10  m 

6, 

=  20.4" 

r  =  268.8  m 

r„  =  100  m 

T 

“-It 

=  150  m 

r,  =  10  km 

z,  =  25  m 

0, 

=  1 35° 

II 

TDPE  can  be  used  to  compute  the  quantity  u,(r  ,z,8'  ,/) 
=  pl{r'.z,6 \t  +  r'/c„),  where 


U,  ( r\z,0 ' )  = 


u,(r',z,6',t)exp(i(ot)dt.  ( 30 ) 


and  p,  is  the  incident  field.  Since  kltL  =  0{  1 ),  we  assume  a 
finite  bandwidth  <  <y  <  w,  and  invert  the  Fourier  trans¬ 
form  in  Eq.  (27)  to  obtain 

(/>vC'=  f  U.ir^e,  +  ir)D,{8jO) 

<  Gm  (x,x„)exp[  —  uo(t  —  r,/c„)  ~\d(o.  (31) 

We  substitute  Eq.  (30)  into  Eq.  (31)  and  interchange  the 
order  of  integration  to  obtain 


.  e  =  -Lf( ft,  -  d—±F)  -  ±f( ft,  - 

d  \  c,,  )  d ,  V  c,.  J 


(P 


d\  =  r  +  (z±z,.)\ 


F(  Oj) 


K(0,t) 


J 


K(0,t  —  t '  )u,(rl,zlpd,  +  ttJ  ')dt ’, 


1  f" 

rAD‘ 


((?,0)exp(  —  i<ot)d(o. 


(32) 

(33) 

(34) 

(35) 


The  integral  operator  in  Eq.  (34)  transforms  incident  wave¬ 
forms  into  scattered  waveforms.  IfEq.  (28)  is  valid,  (p, )'"" 
may  be  used  to  initialize  the  TDPE  at  r  =  r„  <  r( . 

To  illustrate  the  time-domain  solution,  we  consider  the 
Hanning  weighted  sinusoidal  source  function21* 


At)  = 


( 1  —  cos  i<yr)sin  cot, 

0, 


0  <t<T, 
otherwise, 


(36) 


Range  (m)  =  10000.00 


ion  I—  .  . . 1 

BOO  600  400  200  0 

Relative  time  (ms) 


where  co  =  50rrs  1  and  T  =  Stt/oj,  and  an  example  for 
which  data  appear  in  Table  III.  The  incident  field p ,  comput¬ 
ed  with  the  TDPE  appears  in  Fig.  6.  The  ocean  bottom  is 
marked  with  a  horizontal  line.  The  solid  contours  corre¬ 
spond  to  p,  >0;  the  dashed  contours  correspond  to  p,  <0. 
Although  the  time-dependent  wave  equation  does  not  sepa¬ 
rate  in  depth  since  c  depends  on  z,  normal  mode  structure  is 
visible  in  Fig.  6.  The  first  mode  is  present  near  i  —  150  ms, 
the  second  mode  near  /  =  250  ms,  and  the  third  mode  near 
t  =  400  ms.  We  used  the  T-matrix  method  to  calculate  Dr 
The  incident  waveform  u,{r,.z(),0,  +  ~,t)  and  the  scattered 
waveform  F(0,,t)  appear  in  Fig.  7.  Equation  (32)  was  used 
to  calculate  (ps );“".  which  appears  in  Fig.  8. 

The  TDPE  was  applied  in  the  specular  direction  to  ob¬ 
tain  the  scattered  field  ( p ,  ),„„,  which  appears  in  Fig.  9.  Nor¬ 
mal  mode  structure  is  also  visible  in  (pv  )„ul .  However,  the 
second  mode  is  not  excited  because  it  vanishes  near  the  cen¬ 
ter  of  T.  Furthermore,  the  arrival  times  of  the  modes  in  the 
scattered  field  are  qualitatively  different  from  the  arrival 
times  of  the  modes  in  the  incident  field  because  the  scattering 
process  redistributes  the  energy  among  the  modes.'  In  con¬ 
trast  to  the  incident  field,  the  first  mode  arrives  twice,  near 
/  =  200  ms  and  near  t  =  400  ms.  The  third-mode  arrivals 
behave  similarly.  Furthermore,  the  arrival  times  of  the  first 
and  third  modes  overlap. 

V.  CONCLUSIONS 

A  matched  asymptotic  solution  has  been  derived  for  un¬ 
derwater  acoustics  problems  involving  both  propagation 
and  scattering.  The  propagation  and  scattering  aspects  of  the 


FIG  6  The  incident  pulse.  Patterns  of  the  first  three  normal  modes  arc 

clearly  visible  FIG.  7.  The  waveforms  (a)  incident  and  (b)  scattered 
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Range  (m)  =  100.00 


U00  600  400  200  0 

Relative  time  (ms) 


FIG.  8.  The  half-space  scattered  field  ( p ,  )•""  at  r  -  100  m  in  the  specular 
direction  0,  —  45°. 


problem  are  decoupled.  The  model  is  valid  for  scatterers  of 
arbitrary  geometry  and  composition  in  three-dimensional 
environments.  Since  the  small-angle  limit  has  been  used,  the 
model  is  efficient  and  accurate.  In  its  most  general  form,  the 
model  requires  the  solution  of  the  half-space  scattering  prob¬ 
lem  and  the  3DPE.  If  the  scatterer  is  far  from  the  ocean 
surface,  however,  the  half-space  scattering  problem  can  be 
replaced  with  the  free-space  scattering  problem.  If  the  scat¬ 
terer  is  far  from  the  ocean  bottom,  the  solution  of  the  3DPE 
is  not  required  because  the  inner  region  extends  to  a  relative¬ 
ly  large  range.  For  low  frequencies,  the  incident  field  can  be 
approximated  by  a  plane  wave,  and  the  scattered  field  can  be 
approximated  by  the  field  due  to  an  azimuthally  specular 
point  source.  These  approximations  make  low-frequency 
pulse  scattering  problems  practical. 
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